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In a recent letter, Fisher et al. reported the phenomenon of zero-lag long range isochronous 
synchronization via dynamical relaying in systems with delay [Phys. Rev. Lett. 97, 123902 (2006)]. 
They reported that when one has two coupled systems A and C, with delay between them, then the 
introduction of a third element B between A and C will allow them to synchronize even in regions of 
the parameter space where this was not possible without the presence of B. Here we study in detail 
the phenomenon and verify that in all the cases studied (including the ones reported by Fisher et 
al.) this occurs due to the tendency of A and B and B and C to be in antiphase synchronization and 
if A is in antiphase with B and B is in antiphase with C, it will imply that A and C are inphase. We 
show this in coupled quadratic maps, Kuramoto and Rossler oscillators. We also report that there 
is a simpler configuration where the same phenomenon occurs and has nearly the same features of 
the cases studied by Fisher et al. 



INTRODUCTION 

Synchronization phenomenon is widely observed in na- 
ture and in artificial systems. The first description of 
synchronization is believed to have been made by Huy- 
gens. He observed that two pendulum clocks suspended 
in the same wooden beam tend to synchronize in oppo- 
site swingsQ. In more recent studies, synchronization 
has been observed in fireflies, pacemaker cells, networks 
of neurons, etc 01 . In brain activity, it has been found 
that cognitive acts, such as face perception, is associ- 
ated with zero-lag (i-e., isochronal) synchronization of 
gamma oscillations [S]. Other studies have shown neu- 
ral synchronization in widely separated areas of cortical 
regions 0,01 . 

In a recent publication. Fisher et al.\B\ reported the 
phenomenon of zero-lag long range synchronization via 
dynamical relaying. That is, when one has two coupled 
systems, with delay between them, introducing a third 
one will allow synchronization between the two outer el- 
ements in regions of the parameter space where this was 
impossible without the middle element. In their letter 
they showed that this occurs for coupled lasers (in ex- 
perimental and numerical studies), coupled neurons and 
stated that they observed the same phenomenon in a 
large variety of dynamical systems (such as excitable sys- 
tems, oscillators and maps). 

The aim of this paper is to provide an explanation for 
this phenomenon. We find that this occurs because two 
identical systems lose synchronization as the delay be- 
tween them is increased in such a way that they get in 
antiphase synchronization. That is, in the case of oscilla- 
tors their phase separation is tt. When a third element is 
added between those elements, then the middle element 
will be in antiphase with the outer elements and therefore 
the outer elements will be inphase. Here we call the at- 
tention to the fact that it is not any finite delay that will 
cause loss of synchronization between two coupled sys- 
tems, as we will see below. For loss of synchronization 
to occur the delay has to be in a given range of values 



and therefore this phenomenon is only important in that 
region. 

In their letter. Fisher et al. studied bidirectional cou- 
pling between the middle and outer elements. Here we 
show that nearly identical results can be obtained in a 
simpler configuration, that is, the subsystem coupling as 
studied in [7|. We demonstrate our results in quadratic 
mapsfl], Kuramoto0 and coupled Rossler [lol| oscillators. 
In the case of Kuramoto oscillators, we show analytically 
that when all the oscillators are identical, synchroniza- 
tion between the two outer elements always occur, inde- 
pendently of the delay and coupling values. Although the 
synchronization properties do not depend on if the cou- 
pling is bidirectional or of the subsystem type, in both 
cases it is necessary that the parameters of the outer el- 
ements be the same. 

The paper is organized as follows: In the next section 
we study zero-lag synchronization via dynamical relay- 
ing in a discrete time system, that is, the quadratic map. 
Then we turn our attention to this phenomenon in contin- 
uous time systems. In particular, we study the Kuramoto 
and Rossler oscillators. The last section is dedicated to 
discussion. 



QUADRATIC MAPS 

Here we study a system of coupled quadratic maps. 
For two coupled quadratic maps with delay of two time 
steps, we have 



nt+i 



(1 - k)f{x[) + kf{x\-^) 
= {l-k)f{x\) + kf{x\-^), 



(1) 



where / = 1 



^ a G [0, 2], with / mapping the inter- 



val [-1,1] into itself. The case when the coupling strength, 
fc, is zero has been extensively studied in the literature 
and one knows that the logistic map presents period dou- 
bling bifurcations, leading to a chaotic behavior as a is 
increased^. In the coupling scheme studied in ^ our 



2 



system would be 



k)f{A) 

k)f{4) + k[f{x\-') + f{4-')]/2. (2) 



(1 
(1 



Note that in Eq. [T] the delay is of two time steps and 
in Eq. [2] the delay is of one time step, since, as done by 
Fisher et al, one assumes that the delay between ele- 
ments 1 and 2, as well as between 2 and 3 is half of the 
delay between 1 and 3. 

For the system of coupled quadratic maps considered 
above, we noticed that more than one basin of attrac- 
tion is present. Therefore, studying the region of syn- 
chronization for a given initial condition, may not give a 
full picture of the synchronization region. We observed 
that, as in the case without delay 0, the transverse Lia- 
punov exponent correctly determines the region where 
synchronization is possible between two elements of a 
given system. The region of synchronization found in this 
way does not depend on the initial condition used, since 
by construction the transverse Liapunov exponent mea- 
sures the separation (or shrinkage) of trajectories that 
are started nearby to each other. We have adapted the 
method of Benettin et al. ll| to study the tranverse Li- 
apunov exponent of the coupled quadratic maps. If we 
are studying the synchronization between element a and 
b, after the transient dies out, we evolve the orbit of b 
by making it slightly different from the orbit of element 
a. Then we verify how the orbits of a and b differ after 
one iteration. The perturbation is renormalized in the di- 
rection of maximum growth and the process is repeated 
many times. The transverse Liapunov exponent is given 
by the average logarithm (in this paper we use base 2) of 
the growth (or shrinkage) of the perturbation along the 
orbit. In Figs.[T](a) and (b) we show the region where the 
tranverse Liapunov is non negative (shaded area) for, re- 
spectively, elements 1 and 2 in Eq.[T]and elements 1 and 3 
in Eq. [51 We see that the region where synchronization is 
possible (i.e., when the transverse Liapunov exponent is 
negative) is greatly increased with the scheme given by 
Eq. [21 Although the study of synchronization between 
1 and 2 in Eq. [2l is not the main goal of this work, we 
show in Fig. [Ijc) for the sake of comparison, the region 
of non negative transverse Liapunov exponent for their 
synchronization orbit. We see that synchronization be- 
tween 1 and 2 in Eq. [2l occurs in a smaller region of the 
parameter space, when we compare with synchronization 
between 1 and 2 in Eq. [H and between 1 and 3 in Eq. [21 
This shows that 1 and 3 can be synchronized even when 
1 and 2 cannot. 

Our next step is to compare these results with the syn- 
chronization properties of 1 and 3 when the system is 
coupled via a subsystem coupling, as studied in Q for 
systems without delay in the context of chaotic synchro- 
nization. In this situation, the equations that govern the 
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FIG. 1: Region of nonnegative transverse Liapunov exponent 
for synchronization (shaded) between (a) 1 and 2 in Eq.[T] (b) 
1 and 3 in Eq.[2l (c) 1 and 2 in Eq.[2] (d) 1 and 3 in Eq.|3 (e) 
1 and 3 for subsystem coupling without delay and (f) 1 and 2 
for subsystem coupling without delay. The initial conditions 
used were for Eq. [H xj"^ = 0.3, x^^ = 0.1, a;]"^ = 0.25, 
= 0.05, x°i = 0.2 and xl = 0.0. For Eqs. [2]and Eqs. [3 



we used: a;, = 0.25, 



= 0.05, x^^ 



0.55, xl = 0.2, 



x% = 0.0 and xjj = 0.5. For the subsystem coupling without 
delay we used: xl = 0.2, X2 = 0.0, x'i = 0.5. A transient 
of 20000 iterations was discard and the Liapunov exponent 
was considered non negative if it was greater than 10~^ after 
50000 iterations. 



dynamics of the coupled system are 



= {l-k)f{x\) + kf{x\-^) 

= (1 - fc)/(4) + k,f{x\-') 

= (1 - k)f{x\) + k,f{xl-'). 



(3) 



We show in Fig. [T](d) the region of non negative trans- 
verse Liapunov exponent (shaded area) for the orbits of 
elements 1 and 3. When compared with Fig. [TJb) one 
sees that the two shaded regions are nearly identical to 
each other, showing that there is basically no difference 
of the parameter region where synchronization can occur 
between 1 and 3 in one coupling versus the other. We 
have noticed only very minor differences in the parameter 
region of synchronization between 1 and 2 in Eq. [31 when 
compared with Fig. [Ijc). For the sake of completeness, 
we also show in Figs. [He) and (f ) the region of non neg- 
ative Liapunov exponent between 1 and 3 and 1 and 2, 



3 



1 

-1 
-2 

-4 
-5 
-6 
-7 



1 1 1 1 


1 1 1 1 




// 


' //\ 


! 


/ / ' 

/ 


- 


\ 




- (a) 
1 1 1 1 


1 1 1 1 



0.20.40.60.8 1 1.21.41.61.8 2 
a 




t 



FIG. 2: (a) The transverse Liapunov exponent for the syn- 
chronization orbit of 1 and 2 (solid) in the coupled quadratic 
maps for the subsystem coupling. The dashed lines are the 
transverse Liapunov exponent between 1 and 3 for the sub- 
system and bidirectional coupling, (b) Time evolution of xi 
(solid) and X2 (dashed) for parameter values where synchro- 
nization between 1 and 3 occurs. The parameter values are in 
(a) k = 0.2 and in (b) k = 0.2, a = 1.4. The initial conditions 
are the same as in Fig. [1] 



respectively, for a system without delay (for subsystem 
coupling). When compared with the case with delay, 
we see that delay helps in the synchronization of 1 and 
3 (since the transverse Liapunov is negative in a wider 
area) and it is detrimental for synchronization between 1 
and 2. 

We have studied the magnitude of the transverse Li- 
apunov to find out which coupling, bidirectional or sub- 
system, results in faster synchronization. The more neg- 
ative the transverse Liapunov is, the faster is the syn- 
chronization. As shown in Fig. [H^a) (dashed lines) the 
Liapunov exponent is essentially the same independently 
of the coupling type. In Fig. [Ub) we show the tempo- 
ral evolution (for the subsystem coupling) of the orbit 
of element 1 and the orbit of element 2 in the region 
where 1 synchronizes with 3. We see that they are in 
antiphase synchronization, i.e., moving in opposite direc- 
tions. Thus, 1 is in antiphase with 2, which in turn is in 
antiphase with 3 and consequently 1 is inphase with 3. 
So, this is the mechanism by which the relaying element, 
2, helps in the synchronization our the outer elements, 1 
and 3. 



KURAMOTO OSCILLATORS 



The equations for the subsystem coupling in coupled 
Kuramoto oscillators are given by 



Pi — uj — ksin{(pi — yj2) 
i>2 = uj' ~k sin(02 - 4>i) 

/csin(03 - 02): 



(4) 



UJ 



where we have used the notation <jf to denote (j){t — r) 
and 4) to denote 4'{t). We see that oscillators 1 and 2 are 
uncoupled from oscillator 3 and in this way the analysis 
of two coupled oscillators is helpful in the understand- 
ing of the system as a whole. A system of two coupled 
Kuramoto oscillators with delay was studied by Schuster 
and Wagner and is governed by 



bi ~ — /c sin(0i — 02) 
t>2 — Lo' — k sin((/)2 — (pl) 



(5) 



For the sake of completeness, here we reproduce the main 
results of • They assumed that the synchronizing so- 
lution for the two coupled oscillators is given by 



01,2 = V' ± a/2, 



(6) 



where ip is the common dependent phase and a is a con- 
stant phase shift. Inserting this into Eq. [5]one gets 



cos(-!/' — ?/'^) = (cj — a;')/(2fcsin(Q!)) = const, 



(7) 



which is fulfilled ii ip = fit. Thus, the synchronized solu- 
tion is given by 0i,2 = Qt ± a/2. Inserting this solution 
into Eq. [5] we get that the synchronizing frequency is 
given by the zeros of 



f{fl) =tU-f7-/ctan(r2T)Vcos2(rjT) -fc2/A;2, (8) 

where kc = {uj — uj'^jl, uj = (uj\ + uj2)/1. It is found 
that multiple solutions are possible for when delay is 
present 12] . This is in contrast with the case without 
delay. The phase shift that belongs to the synchronized 
solution is 

a = arcsin{kc/ k cosiyiT)) , if cos(r2T) > 0, 

a = TT — arcsm(fcc/A: cos(57t)), otherwise. (9) 

li UJ = uj' the above equations will reduce to 

f{n) ^uj -n-kfi\n{VLT), (10) 

and 



As in the case of discrete time system, we have found 
that the synchronization properties of the subsystem cou- 
pling are nearly identical to the synchronization proper- 
ties of the bidirectional coupling. Because the subsystem 
configuration is simpler and allow us to perform some 
analytical studies in the Kuramoto system we will con- 
centrate from now on mostly on the subsystem case. 



a = 0, if cos(riT) > 0, 

a = -K. otherwise. (11) 

In Fig. [3ja) we show the region where synchronization 
between oscillator 1 and 2 occurs (white) and when it 
does not occur for a system with uj = uj' = \. The phase 
separation in this case is tt as shown above. The region 
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FIG. 3: Region where isochronal synchronization occurs 
(white) between oscillator 1 and 2 for the Kuramoto oscil- 
lators and when their phase separation is vr (shaded) for 
(a) subsystem coupling and (b) bidirectional coupling, with 
uj = u' — 1. Region where isochronal synchronization occurs 
(white) between oscillator 1 and 3 for the subsystem coupling 
with (c) Lo = 1.2, uj' = 1, (d) LO ^ 1, uj' ^ 1.2. The initial 
conditions used were (pi = 0.1, (j)2 ~ 0.12, 03 = 0.14 and 
(/>i,2,3(-r) = <^i,2.3(t = 0) +0.001r. 

of non synchronization is in the shape of " tongues" with 
its center being around n7r/2, where n is an odd number. 

We have found that when uj = uj', oscillator 1 and 3 
always synchronize, independently of the delay and cou- 
phng values. We can show this in the following way: give 
a small perturbation A, where A is positively defined, to 
the orbit of oscillator 1, so that the initial condition for 
oscillator 3 is given by (/)3 = 0i + A. Now subtracting 
the first from the third lines of Eq. 3] we have 

(p3-<j)i^ k[sm{(f>2 - 03) - sin(02 - M], (12) 

Using A as defined above and the Taylor series expansion 
we get from Eq. [7] after a straight forward algebra, 

A = -fcAcos(02 - 01) = -fcAcos(17T - a), (13) 

We have seen that if cos(f2r) > 0, then a = and 
if cos(r2T < 0), then a = tt. Both cases result in 
A/A = — fc|cos(r2T)| < 0. Consequently, the separa- 
tion A between the orbits of oscillator 1 and 3 decreases 
as time evolves, implying that they tend to synchronize 
for any value of the parameter space (with A; > 0). We 
also notice that there is region of the parameter space 
where the phase differences (say, or tt in the case of 
identical oscillators) do not change if the delay is varied 
in that region. Consequently, we expect that even if the 
delay between 1 and 2 is not exactly the same as between 



2 and 3, zero-lag long range isochronal synchronization 
between 1 and 3 would still be possible. 

For the sake of comparison, we show in Fig. the 
region of synchronization between 1 and 2 (white) for the 
bidirectional coupling, that is, 

01 — — fc sin(0i — 02) 

02 = c^'-fc[sin(02-0I)+sin(02-05)]/2 (14) 

03 = w - fcsin(03 - 02). 

We see that Figs. [3] (a) and (b) are almost identical, con- 
firming once more that there is basically no distinction 
between the synchronization properties of the subsystem 
and bidirectional couplings. 

We have also considered the case in which lu uj'. 
This is displayed in Figs. [3] (c) and (d). When to ^ uj' 
there is no isochronal synchronization between 1 and 2, 
since there is a constant phase difference between the 
oscillators [l^. However, synchronization between 1 and 

3 is possible in a large area of the parameter space. The 
direction of the "tongues" where synchronization does 
not occur depend on having lo > lo' ot the other way 
around. However, the area of non synchronization varies 
little from one case to the other. 

We have looked for even simpler configurations where 
the zero-lag long range isochronal synchronization would 
occur, such as removing one of the links between element 
1 and 2 in the subsystem coupling, but we found that 
such a kind synchronization is not possible in those sim- 
pler configurations. 

ROSSLER OSCILLATORS 

We have chosen Rossler chaotic oscillators to show that 
also in this chaotic system, governed by continuous time 
equations, the phenomenon reported in Q is due to the 
fact that inphase synchronization coupled with another 
inphase synchronization results in isochronal synchro- 
nization. The equation of Rossler oscillator with sub- 
system coupling are given by 

Xi = -Loyi ~ Zi- k{xi - Xj) 

iji = LOXi + ayi - k{yi - yj) (15) 

Zi = b + z{xt - c) - k{zi - zj), 

where j = 1, 2, 3, and when i = l,j = 2, i = 2,j = l and 
i — 3, j = 2. In this paper we have used, cj = 1, a = 0.15, 
b = 0.4 and c = 8.5. 

It is shown in Figs. [H^a) and (b) the region where syn- 
chronization happens (white) between oscillators 1 and 2 
and between oscillators 1 and 3, respectively. We see once 
more that the presence of a middle element increases the 
region where synchronization is possible. In Fig. [DJc) 
we show the time evolution of yi{t) (solid) and y2{t) 
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FIG. 4: Region where synchronization occurs (white) and 
when it does not (shaded) between (a) 1 and 2 and (b) 1 and 
3 for the Rossler coupled oscillators with subsystem coupling, 
(c) Temporal evolution of yi{t) (solid) and y2{t) (dashed) 
when T — 2.5 and k — 0.2. The parameter values were 
LU = 1, a — 0.15, b = 0.4 and c = 8.5, and the initial con- 
ditions where xi = 0.1, yi — 0.12, zi — 0.14, X2 = 0.1, 
3/2 = 0.22, Z2 = 0.24, X3 = 0.1, ys = 0.12, zi = 0.14 and 
x,y,z(—T) = x,y,z{t — 0) + O.OOlr. (c) Temporal evolution 
of xi{t) (solid) and X2(t) (dashed) when r — 2.5 and k = 0.2. 



(dashed) for a region where oscillator 1 and 3 synchro- 
nize, that is, T — 2.5 and k — 0.2. Once again we see 
the antiphase synchronization observed in the previous 
sections. 



DISCUSSION 

By studying the results presented in the paper by 
Fisher et ai., we see the origin of for zero-lag long range 
isochronal synchronization they observed is the same we 



have found. If we take a look on Figs. 2(E), 2(F), 3(E) 
and 3(F) of @ we see that they shifted the time series 
of the middle element by the delay value in order to fa- 
cilitate the comparison with the time series of the outer 
elements. When they did this shift, the time series all 
the elements nearly matched each other. But it turns 
out that, apparently not realized by them, the value Tc is 
the result of this antiphase mechanism mentioned in the 
previous paragraphs. In our paper we have clarified this 
point, and showed that there is a simpler configuration 
(the subsystem coupling) where nearly identical proper- 
ties are observed when compared with the bidirectional 
coupling. We have also been able to perform analyti- 
cal studies and show that for the Kuramoto system with 
identical oscillators, synchronization between the outer 
elements always occurs, independently of the parameter 
values. 
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